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ABSTRACT 



We have performed C^^O (J=l-0) mapping observations of a 20' x 20' area 
of the OMC-1 region in the Orion A cloud. We identified 65 C^^O cores, which 
have mean radius, velocity width in FWHM, and LTE mass of 0.18±0.03 pc, 
0.40±0.15 km s~^, and 7.2±4.5 Mq, respectively. All the cores are most hkely to 
be gravitationally bound by considering the uncertainty in the C^^O abundance. 
We derived a C^^O core mass function, which shows a power-law-like behavior 
above 5 Mq. The best-fit power-law index of —2.3 ± 0.3 is consistent with those 
of the dense core mass functions and the stellar initial mass function (IMF) 
previously derived in the OMC-1 region. This agreement strongly suggests that 
the power-law form of the IMF has been already determined at the density of 
~ 10^ cm~^, traced by the C^^O ( J=l-0) line. Consequently, wc propose that the 
origin of the IMF should be searched in tenuous cloud structures with densities 
of less than 10^ cm~^. 

Subject headings: ISM: clouds — ISM: individual(Orion A) — stars: formation 
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Introduction 



One of the most important observational features of the stellar initial mass function 
(IMF) is its power-law-like nature above 1 Mq, as dN/dM oc M~^. In the solar 



1955 



Kroupa 



neigh borhood, the power-law index 7 seems to be greater than 2 (iSalpeterl 
2OOII ). which characterizes the statistical properties of stars. In particular, both the total 
number and mass of stars are dominated by those of low- mass stars of ~ 1 Mq. It is 
natural to consider that the origin of the IMF shape is related to the mass distribution of 
its natal gas in molecular clouds. Many authors have investigated dense gas (10^~^ cm~^) in 
molecular clouds by using (sub) millimeter dust continuum emission and/or molecular line 



emiss ions having 



ligh c r itical densities such as the H C0"^( J=l-0) and ^f:^H 



(e.g. 


Motte et al. 


1998; 


Reid & Wilson 


2006 


2007; 


Walsh et al. 


2007 


Enoch et al. 


2008 


)• 



Nutter fc Ward- Thompson 



f J=l-0) lines 



2007 



Ikeda et al. 



typical sizes of 0.05 - 0.1 pc and masses of 1 - 10 Mq in nearby (< 500 pc) star forming 
regions such as Orion, Ophiuchus, Perseus, and Serpens. The molecular line studies showed 
that the cores are gravitationally bound and are likely to produce stars. Moreover, they 
found that the core mass functions (CMFs) derived by using the dense gas tracers, referred 
to as DCMFs hereafter, seem to have power-law-like behaviors in hi gh-mass parts , whos e 



Kramer et al. 



fll998h 



7 values are very similar to that of the IMF. One exception is that 
found a significatly smaller 7 value of 1.7 for the S140 and M17SW regions by using 
the C^^0(J=2-1) line with a high critical density comparable to that of II^^CO^(l-O). 
Considering that the power-law form of the IMF has been already determined at the 
formation stage of the cores with the densities of 10"^"^ cm~^, it is likely that more tenuous 
structures of molecular clouds have a key to understanding the origin of the power-law 
nature of the IMF. 



It has been suggested that the mass functions in the tenuous gas structures are different 



-4- 



from the DCMFs. iKramer et al.l (119981 ) carried out a systematic study of the mass functions 
in the tenuous gas structures of 10^ cm~^ or less by using the ^^C0(2-l), ^^CO(1-0; 2-1), 
and C^^O(l-O) hues in various molecular clouds. They showed that their mass functions 
seem to have a common power-law for m, and the 7 value of 1.7 ±0.1 is significantly smaller 



than those of the DCMFs and the IMF. 



Heithausen et al. 



(119981 ) derived ^'CO(J=1-0, 2-1) 



mass fun ctions in the MCL D 123.5+24.9 and Polaris Flare regions and found that 7 



1.8±0.1. 



Wong et al. 



((20081) found 7 = 1.7 in the 0^^0(7=1-0) mass function of ROW 106. 
These facts mean that the power-law index of the IMF may be originated in the formation 
process of the dense gas of 10^~^ cm~'^ from the tenuous gas of 10^"'^ cm~^. However, one 
should be careful in comparing the 7 value of the tenuous gas mass function to those of 
the DOMFs and the IMF. This is because the tenuous gas mass functions described above 
were derived by the spatial resolutions larger than 0.1 pc, which cannot resolve star-forming 
cores, and/or by using optically-thick tracers such as ^^00 and ^'^00. To fairly compare 
the tenuous gas mass function with the DOMFs and the IMFs, one should achieve higher 
spatial resolutions than 0.1 pc and use optically-thin tracers. 

In this paper, we present a OMF derived by 0^*^0(7=1-0) mapping observations of 
the OMO-1 region, which is located at the center of the Orion A Oiant Molecular Oloud. 
The aim of this study is to examine whether or not the common power-law form between 
the DOMFs and the IMF has been already determined in the tenuous gas of < 10^ cm~^ 



molecular line emission is suitable for deriving the OMF (e.g.. 



Tachihara et al. 


2002 


) 


1 

~^ ((Ungerechts et al. 


1997) 



and is typically optically thin (see §3.2p . The OMO-1 region is one of the best regions to 



investigate th e OMF, because t 



l e IMF of the associa ted Orion Nebula Oluster (ONO) has 



1997 



Muench et al 



2OO2I ) ■ and because the power- 



been derived ( [Hillenbrand 
the H^300+ DOMF is shown to be quite similar to that of the ONO IMF by 



aw form of 



Ikeda et al 
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( 120071 ) ■ Furthermore, at the distance to the Orion A cloud of 480 pc (iGenzel et al.lll98ll ) 
we can easily resolve the cores with radii of ~0.1 pc. As described in §2], the mapping 
observations have been done with the effective spatial resolutio n of 26^^4 (= 0.06 pc), which 



is high enough to resolve the dense cores in the OMC-1 region (jlkeda et al. 



20071). 



Observations 



The mapping observations were carried out in the period from January to February 
2008 by using the Nobeyama 45 m radio telescope. The mapping covered a central 20' x 20' 
area of the OMC-1 region. The position of the Orion-KL object (5^35™14^2, -5°22'22"; 
J2000) was taken to be the center of the mapping region. Note that the regi on was selected 



to be the same as th at in the H^^CO"*" D CMF study by 



On-The-Fly method (iSawada et al. 



Ikeda et al. 



(120071 ). Using the 



20081 ). we swept the mapping region by raster scan 



with a scan speed of the telescope of 52 arcsec s ^. To reduce scanning effects, we scanned 
in both the RA a nd Dec directions. At the frequency of the C^^O(J=1-0) emission 



(109.782182 GHz; 



Ungerechts et al. 



19971 ). the half power beam width, A^hpbw, and main 



beam efficiency, rj, of the telescope were 14" and 0.4, respectively. At the front end, we 
used the 25-BEam Array Receiver System (BEARS) i n double-sideband r node, which has 



5x5 beams separated by 41". 1 in the plane of the sky (ISunada et al. 



2000 



Yamaguchi et al. 



20001 ) . The 25 beams have beam-to-beam variations of about 10% in both beam efficiency 
and sideband ratio. To correct for the beam-to-beam gain variations, we calibrated the 
intensity scale of each beam using a 100 GHz SIS receiver (SlOO) with a single-sideband 



filter. At the back end, we used 25 sets of 1024 channel autocorre 



ators (AGs), which 



20001 ). Since the data 



have a velocity resolution of 0.104 km s~^ at 110 GHz (jSorai et al.l 
dumping time of the AGs was 0.1 s, the spatial data sampling interval on the sky plane 
was 5.2". The interval corresponds to 0.35A^hpbw5 and satisfies the Nyquist theorem. The 
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telescope pointing was checked every 1.5 hours by observing the SiO (f=l, J=l-0; 43.122 
GHz) maser source Orion KL. The pointing accuracy was better than 3". To construct an 
LSR data cube, we used a Gaussian function as a gridding convolution function (GCF) 
to integrate the spectra which were taken with a very high spatial sampling rate of 5.2". 
We adopted 22". 5 as the size of the GCF in full width at half maximum. The resultant 
effective spatial resolution of the cube, A6'efr, becomes 26". 4. During the observations, the 
system noise temperature ranged from 295 to 429 K with a mean value of 351 K. Therefore, 
we have an RMS noise level of the data cube of 0.18 K in T^. 



Results 



3.1. 0^^0(7=1-0) Total Map 

Figure [1] shows the total integrated intensity map of the C^^O(J=1-0) emission. The 
map shows that the OMC-1 region is elongated along the north- south direction. T his 
elongated structure is a part of the Integral Shaped Filament (ISF; 



Bally et al 



1987), which 



Dutrey et al. 



i s a g lobal structure of the northern part of the Orion A cloud. Previously, 
(Il99ll ) obtained the 0^*^0(7=1-0) map of the ISF with a spatial resolution 6 times coarser 
than ours. Although the overall feature of our map is consistent with theirs, our map 
reveals a number of 0.1 pc-scale sub-structures owing to our higher spatial r esolution. Su ch 
a clu mpy structure is also seen in the H^^CO"^ map of the OMC-1 region (jlkeda et al 



20071 ). which traces the dense gas of ~ 10^ cm ^, one order of magnitude higher than ~ 10'^ 



cm ^ traced by the C^^O line. 

One distinct difference between the C^^O and H^^CO^ maps is that the C^'^0 map 
seems less clumpy than the H^^CO^ one, though the overall features of the two maps are 
similar to each other. In the H^^CO^ map, the ISF seems to mainly consist of the 0.1 
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pc-scale cores and the large-scale diffuse component seems minor. On the contrary, the 
diffuse component dominates the C^^O map, and the clumpy structure seems less distinct 
than that in the H^'^CO^ map. In addition, the diffuse component of the C^^O map appears 
to be more extended than that of the H^^CO^ one. Actually, the total mass of the C^^O 
emission above the 2 a level is estimated to be 1500 Mq (see §3.21 for the mass estimate) 
and is 1.5 times larger than that of the H^^CO+ emission in the same region (~ 1000 Mq] 



Ikeda et al. 



20071 ). These facts support that the C^^O emission traces more tenuous gas of 



the cloud than the H^3co+ emission. 



3.2. Identification of C^^O Cores and Estimation of the Core Properties 



Following the Hi3C0+ DCM F study in the OMC- 1 region fllkeda et al. 



applied the clumpf ind algorithm (iWilliams et al. 



2007|), we 



1994J ) to the C O three-dimens i onal 



(a-5-fLSR) data. For the input parameter for the algorithm, AT in 



Williams et al. 



fll994h . 



we adopted 0.36 K, i.e., the 2 a noise level of the cube data. Since the algorithm treats 



the contour levels of nAT with n 
contour level of AT. Williams et al. 



, 2, 3 , ■ ■ ■, the threshold level is equal to the lowest 



AT to extract clump and/or core structures from t 
properties of them, including the CMF. R ecently, 



(1l994l ) showed that the 2 a noise level is optimal for 
le cube data a nd to recover the physical 



wider ranges of AT (from 3 to 20 a) than 



Pineda et al 



Williams et al 



( 120091 ) examined more 



(I1994I ) did and found that the 



mass functions of cloud sub structures specific to the adopted AT have power-law indices 
depending on AT, especially for higher AT (> 5 a), suggesting the hierarchical nature of a 
molecular cloud. In this study, however, we focus our attention on the cores, corresponding 
to the 0.1 pc-scale hierarchical level directly related to star formation. Therefore, we prefer 
to use the clumpf ind algorithm with the optimal AT of the 2 a noise level confirmed by 



Williams et al 



(1199! ^see also Note that we adopted the grid spacing of the cube 
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data o f 26". 4, equal to A^efr, i-e., full-beam sampling. This is because IWilliams et al. 

I-bea m sampling case. We 



( 119941 ) determined the optimal AT of the 2 a le vel for the fu . 



Ikeda et al 



20071 ) to reject ambiguous or 



also followed the additional criteria introduced in 
fake core candidates whose size and velocity width are smaller than the spatial and velocity 
resolutions, respectively. As a result, we identified 65 cores. The total number of the C^^O 
cores is comparable to that of the H^^CO+ cores of 57 in the same region. 

We estimated the radius Rcore, velocity width in FWHM dvcoie, LTE mass Mlte, virial 
mass Mvir, and mean density n o f the cores. The definitions of these parameters are 

the same as those in 



Ikeda et al 



fl2007f ). In this study, A^eff = 26".- 



km s ^. In the mass estimate, we adopted rj = 0.4 and Tgx = 20K (jCesaroni fc Wilson 



and dv, 



0.10 



994h. For the 



( iFrerking et al. 



r acti onal abundance of C O relative to H2, Xqisq, we adopted 1.7 x 10 



-7 



19821 ). In addition, we confirmed that the 0^*^0(7=1-0) emission is optically 
thin. We observed the ^^OO(J=l-0) emission (110.201 GHz) toward 25 positions in the 
OMC-1 region including the most intense O^^O peak at the Orion-S object (a =5'^35™13'^, 
6 =— 5°24'30"). We derived the optical depth of the O^^O emission tqisq by considering the 
terrestrial abundance ratio of [^^0^^0]/[^^0^^0] of 5.5, and found that tcisq <^ 1 at all the 
25 positions. Therefore, assuming the O^^O emission is optically-thin all over the observed 
area, we have 

-1 



Mlte = 5.0x10 



-2 



X, 



CI80 



1.7 X 10 



-7 



rp ^5.27/Te: 



D 



480pc 



26".4 



_r]_ 
0.4 



-1 /sjXi^ 

K km s~^ 



T. 



where SjT^ -Avi is the total integrated intensity of the core. 



The mean values with standard deviation of -Rcore, dvcove, and Mlte are 0.18±0.03 pc, 
0.40±0.15 km s~^, and 7.2±4.5 Mq, respectively. These are consistent with those of the 
H^^CO+ cores of 0.15±0.04 pc, 0.58±0.23 km s"\ and 15±16 Mq. The mean value of n of 
(4.8±1.6) X 10^ cm~'^ is comparable to the critical density of the C^^O(J=1-0) line and is 
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one order of magnitude smaller than that of the H^^CO^ cores of (3.3±1.8)xl0^ cm ^. The 
total mass of the C^^O cores is 468 M©, 30% of the total mass traced by the C^^O emissio n 



see ^3.ip . This fraction is a half of that for the H^^CO^ case of 60 % (llkeda et al. 



2007 



20091 ). supporting that the C^^O emission is dominated by the diffuse component, compared 
with the H^^CO"*" emission. 

All the C^^O cores are likely to be gravitationally bound and have the potential for 
forming stars as well as the H^'^CO^ cores. Figure [2] shows that the C^^O cores, except for 
one core with the smallest mass of 1.3 Mq, are distributed around the line of Myir/MLTE = 1 
as well as the H^'^CO^ cores. The virial ratio, Mvir/MLTE? has a mean value of l.liO.S and 
the maximum value of 2.6. Even if we take the maximum value, the C^^O cores can be 
under virial equilibrium by considering the uncertainty in Xq^sq of a factor 3. On the other 
hand, the minimum-mass core has the smallest virial ratio of 0.03, indicating certainly 
self-gravitating. Although the core is located around the map center of a =5*^35^13^, 
S = — 5°23|19^^ t he core is just on the border of the core identification (see §3.21 and 
Ikeda et al.l 120071 ) and might have large uncertainties in its physical parameters. Even if the 
core would be an unresolved one, it never affects our discussion on the power-law nature in 
the high-mass part of the CMF, as in §11 



4. C^^O CMF and Comparison with the DCMFs and the IMF 

Figure [3] shows the C^^O CMF. The CMF has a turnover at around 5 Mq, and a 
power-law like shape in the high-mass part above the turnover. Above 5 Mq, we applied 
a single power-law function by considering the statistical uncertainties and found that 
the best-fit power-law index 7 is 2.3±0.3. Furthermore, we confirmed that the 7 value is 
insensitive to the input parameter of the clumpf ind, AT. We checked the algorithm by 
changing the step size (i.e., the contour interval) and the threshold level (i.e., the lowest 
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contour level) independently, as shown in Figure HI Note that the parameter ranges were 
limited so that the number of the identified cores could be large enough to construct the 
CMF. Figure H] clearly shows that all the 7 values are consistent with the best-fit value of 
2.3 within the uncertainties. 



Before the comparison with the DCM Fs and the 



Ikeda et al. 



MFlwe correct the C^^O CMF 



(120091 ). The mass density of the 



for the confusion effect using the model by 
low-mass cores of < 5Mq, pcube, in the C^^O cube data is a key parameter in the confusion 
model and is estimated to be 40.0 M© pc^^ km^^ s. The pcube value is consider ably smaller 
than that for the H^^CO+ cores in the OMC-1 region of 99.0 M© pc"^ km~^ s fllkeda et al 



20091 ). This difference probably corresponds to the fact that the C^^O map is less clumpy 
than the H^^CO^ one, as stated in §3.11 The total mass misidentified by the confusion effect 
is estimated to be 123 Mq, and the confusion-corrected CMF is reconstructed by adding 45 
misidentified cores to the low-mass part below 5 Mq. Here we exclude the minimum- mass 
bin of 2 Mq only for the minimum-mass core, because the detection of the minimum-mass 
core could be marginal as stated in §3.21 and because the statistical uncertainty of the 
mass bin is large compared to the other bins below the turnover. Note that the confusion 
effect does not considerably affect the shape of the high-mass part of the CMF, even if we 
include the minimum-mass bin in the correction. In other words, our discussion about the 
power-law index 7 in the high-mass part of the CMF is not considerably changed by the 
confusion effect. 

In Figure [3] we also show the confusion-corrected C^^O CMF. Although the shape 
in the low-mass part of the corrected CMF is considerably changed, the 7 value in the 
high-mass part is consistent with that in the observed CMF within the uncertainties; we 
applied a single power-law function to the corrected CMF and obtained the best-fit 7 value 
of 2.4 ±0.2. 
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We have the 7 value of the corrected C^^O CMF of 2.4±0.2, which is quite consistent 



with those of the DCMFs. The 7 va^ 



effect was measured to be 2.2 ± 0.1 (llkeda et al 



ue of the h 



^■^C0+ DCMF corrected for the con 



2003). 



Nutter fc Ward- Thompson! (120071 ) 



usion 



derived the 850 /im dust continuum DCMF in Orion and found that the best-fit 7 is 



2.2±0.2. Furt 



lermore, our 7 value agrees well with that of the IMF within the uncertainties. 



Muench et al. 



(120021 ) derived the IMF of the Trapezium cluster, which is the center portion 
of the ONC, from the i^'-band luminosity function and found that 7 = 2.2, as shown in 
Figure [3l The agreement between the C^^O CMF and the IMF suggests that, at least in the 
OMC-1 region in the Orion A cloud, the power-law form of the IMF with 7 > 2 has been 
already determined at the formation time of the tenuous structure with density of ~ 10^ 



cm 



In future works, we propose observations using molecular lines having lower critical 
densities, such as ^^CO(J=1-0) and ^^CO(J=1-0), or thermal dust continuum emission. 



2003) 



Telescopes in orbit with spatial resolution of < 1' such as the Akari ([Murakami et al. 
and Herschel satellites would be preferable to seek the origin of the IM F in the tenuous 



structures of molecular /atomic interstellar clouds (e.g.. 



Toth et al. 



2000). 
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35m40s 35m20s 35m0s 5h34m40s 

Right Ascension (J2000) 



Fig. 1. — Total integrated intensity map of the C^^O{J = 1 — 0) emission ( i'lsr = 0.2 — 15.0 
km s"^ ) in the OMC-1 region. The contour intervals are 0.74 K km s~^ (corresponding to 
2a) starting at 0.74 K km s~^. The gray scale bar is shown at the right-hand side of the 
panel. 
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Fig. 2. — Virial ratio (Mvir/MLTE)-iiiass relation of the C O cores (filled 



circles show the same relation of the H^'^CO^ cores in the OMC-1 region (jlkeda et al. 



circles). The ope n 



2003). 
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Fig. 3. — Observed and confusion-corrected C^^O CMFs are shown by the thick-color filled 
and thin-color open circles, respectively. The error bars indicate the statistical uncertainty 
of V^, where is the core number in each bin. The thick- and thin-color solid lines shows 
the best-fit power-law functions in the high-mass part above 5 Mq for the observed and 
corrected CMFs, re spectively (see jSj). The slope in the high-mass part of the IMF of the 



Trapezium cluster (iMuench et al 



2002) is also shown by the dashed line. Note that the 



dN/ dM value of the IMF is arbitrary scaled. 
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2 a threshold level 2 a step size 




1.5 2 2.5 3 3.5 4 4.5 5 1.5 2 2.5 3 3.5 4 4.5 5 5.5 
Step size [a] Threshold level [a] 



Fig. 4. — Number of the identified C^^O cores (thin-color symbols) and power-law index 
7 of the resultant CMF (thick-color symbols) are shown as a function of the step size and 
the threshold level in the left and right panels, respectively. The y-axis on the left- and 
right-hand sides of each panel indicate the 7 values and the core numbers, respectively. The 
vertical error bars with the 7 symbols correspond to the uncertainties in power-law fitting. 
The horizontal dashed line shows the 7 value obtained by the optimal parameters of the 2 
a step size and the 2 a threshold level (see Figure [3]). 



